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RESONANT EQUILIBRIUM CONFIGURATIONS IN 
QUASI-PERIODIC MEDIA: RAM THEORY 

RAFAEL DE LA LEAVE, XIEENG SU, AND LEI ZHANG 


Abstract. We develop an a-posteriori KAM theory for the equilibrium 
equations for quasi-periodic solutions in a quasi-periodic Erenkel-Kontorova 
model when the frequency of the solutions resonates with the frequen¬ 
cies of the substratum. 

The KAM theory we develop is very different both in the methods and 
in the conclusions from the more customary KAM theory for Hamilton¬ 
ian systems or from the KAM theory in quasi-periodic media for so¬ 
lutions with frequencies which are Diophantine with respect to the fre¬ 
quencies of the media. The main difhculty is that we cannot use transfor¬ 
mations (as in the Hamiltonian case) nor Ward identities (as in the case 
of frequencies Diophantine with those of the media). 

The technique we use is to add an extra equation that ensures the 
linearization of the equilibrium equation factorizes. To solve the extra 
equation requires an extra counterterm. We compare this phenomenon 
with other phenomena in KAM theory. It seems that this technique could 
be used in several other problems. 

As a conclusion, we obtain that the perturbation expansions devel¬ 
oped in the previous paper BdlLSZlSIl converge when the potentials are 
in a codimension one manifold in a space of potentials. The method of 
proof also leads to efficient (low storage requirements and low operation 
count) algorithms to compute the quasi-periodic solutions. 

Quasi-periodic Frenkel-Kontorova models, resonant frequencies, equi¬ 
libria, quasi-crystals, Lindstedt series, counterterms, KAM theory 
[2000] 70K43, 37J40, 52C23 


1. Introduction 

The goal of this paper is to develop a KAM theory for the functional 
equation: 

(1) viijj -I- Q) -I- v(i/r - G) - 2v(i/^) + W{(il/, rf) + -f d = 0 
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where > R, Q g R^^\ yS 6 R^ are given, is a variable in and 

we ean think of 77 G T as a parameter. We are to find v : x ^ R as a 

function of (t/r, 77) and d G R as a function of the parameter 77. We will refer 
to O as the “equilibrium equation”. 

The equation (H]) was derived in the paper HdlLSZlSIl as the equation is 
satisfied by hull functions of quasi-periodic equilibria in Frenkel-Kontorova 
models of deposition on quasi-crystals when the frequency of the equilib¬ 
rium solution is resonant with the frequencies of the substratum. The vari¬ 
able 77 is an angle which has the meaning of a transversal phase. 

Roughly, the model describes particles interacting with their neighbors 
and with a substratum which is quasi-periodic. The configuration describ¬ 
ing the state of the system is parameterized by the hull function v. We try 
to place the particles with a frequency (inverse of the density of particles) 
which resonates with the frequencies of the medium. Q represents the in¬ 
trinsic frequencies. Since the medium is resonating with the frequency of 
the configuration, the positions of the particles are parameterized hy d - I 
angles, i.e. they cover densely did - I dimensional torus which is indexed 
by 77. The W represents the forces of the particles with the substratum and 
the A is an external force. 

We refer to the paper HdlLSZlSlI for the discussion of the physical mo¬ 
tivations (there are several physical motivations for the Frenkel-Kontorova 
model) and for a formal analysis. From the strictly logical point of view, this 
paper and HdlLSZlSlI are completely independent even if they are motivated 
by the same physical problem. They also rely on very different techniques. 
To avoid repetition, we refer to HdlLSZlSl for references to the previous 
literature on the problem as well as for physical motivations. 

The main goal of this paper is to develop a KAM theory for the equation 
O, but we will have to add a one dimensional extra parameter to it. 

The main source of the difficulty to implement a Newton method - as 
needed in KAM theory - is that the equation (l3^ (the linearization of the 
equilibrium equation ([T])) is not easily analyzable in a way that leads to 
tame estimates. We will deal with this problem by adding an extra auxiliary 
equation which implies that (l3^ can be solved with tame estimates. The 
addition of an extra equation that allows to solve the linearization is simi¬ 
lar in spirit to the introduction of the reducibility in KAM theory [|Mos67ll . 
Nevertheless, our auxiliary equation is very different from the one in re¬ 
ducibility. 

1.1. The method of adding extra parameters to equations. The main 
observation that allows us to develop a KAM theory is that if we are allowed 
to adjust a one dimensional parameter in the potential, then the linearized 
equilibrium equation admits a very nice structure (it can be factorized into 
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two first order equations and the faetorization allows to develop an iterative 
proeedure whieh is quadratie eonvergent). 

Henee, in Seetion ILLll we will add an extra parameter to the left hand 
side of dU) to obtain the modified equilibrium equation @ sueh that its lin¬ 
earization ean be faetorized. Then, we supplement Q with another equation 
dll) in Seetion 11.1.21 (we eall it the factorization equation) whieh eneodes 
that the linearization of the equilibrium equation ean be solved. 

1.1.1. The modified equilibrium equation. For eaeh fixed q e we will 
look for a funetion v : > R and for numbers T, cr in sueh a way that 

we have 

q) 

(2) =v(i/f -I- Q) -I- vifif - Q) - 2v(i/r) -I- q) + fiv(fi)) + v{ij/)a + T 
= 0 . 

In the rest of the paper, we will eall this modified equation dZ]) the equilib¬ 
rium equation. The equation ^ has a symmetry that makes the solutions 
not unique (this eorresponds to a gauge symmetry related to the ehoiee of 
origin of the phase in the original problem). Henee, to obtain loeal unique¬ 
ness, we supplement d3) with the following normalization: 

(3) I v{ifj) dip = 0. 

We eonsider (f as a funetional that given a funetion v and two numbers 
cr, A produees another funetion given by the seeond line in dl]). And we treat 
the equation dZ]) as searehing for zeros of the funetional S’. 

Sinee v, cr, A all depend on the parameter q, for eonvenienee, we will write 
Vri when we need to emphasize the faet that v depends on q, and similarly 
for all the other funetions. We will also obtain cr, T as a funetion of q (and 
maybe of frequeney Q, but we will not diseuss dependenee on Q in this 
paper). 

Depending on the physieal solutions, we may impose the value of one of 
these variables and determine the others. For example, if we are imposing 
an external foree we may want to fix A or if the material is eonstrained to 
have a eertain density, we may fix D. 

It is important to notiee that, onee we have established the KAM theorem, 
eliminating some variables in terms of the others is just an applieation of the 
finite dimensional implieit funetion theorem. 

Remark 1. Variants of the idea of adding external parameters and then 
setting them to zero, has appeared in many guises. In perturbative expan¬ 
sions in Physics, it is called the method of eounterterms HBSSOl IGalSSH . 
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In differential equations, it is called Cesari’s alternative method (Chapter 
IX of HHalSOI , [ICH82B ). Closer to us, in KAM theory, it was introduced in 
llMos67ll . It was realized in lldlL861 [BHS961 [^e92l ISev99ll that it provided 
a good way to deal with degenerate problems. A very systematic treat¬ 
ment of dependence on parameters (including parameters taking values in 
nowhere dense sets) appears in [IVan02ll . 

Remark 2. The exact form of the counterterm added is not that important. 
We could have put other counterterms o'F(v)for almost any function F ^ 0. 

The way of thinking geometrically of the counterterms is that there is 
codimension 1 set of potentials for which the solutions move differentiably 
with respect to parameters. The counterterm is a projection that moves to 
keep the problem in this manifold. We could have taken any other family of 
corrections to the codimension one manifold where the solutions are found. 
We refer to Section\^or a precise formulation and more details. 

Remark 3. Even if we could consider Q as afunctional equation for each 
value ofq, we will show in the next section that the symmetries of the equa¬ 
tion involve mixing the f and q dependence. Relatedly, we note that the 
equation ^ for a fixed value ofq does not have a variational principle. 

It will be important to mention that, because jS has components both in 
the and the q directions, the equation Q cannot be considered just as a 
parameterized version of the equations considered in IISdlL12bll . 

1.1.2. The factorization equation. For the KAM treatment of the equilib¬ 
rium equation Q, we will find it useful to supplement ([21) with another 
equation whieh we eall the faetorization equation 

(4) ^[v,cr,c](i/r,77) = {-c{fi)+2-diiW{{il/,q)+l3vfifi))-o-]c{fi+Q.)-\ = 0. 

This should be eonsidered as a funetional equation for funetions v, c : —s 

R and number cr when q is fixed. 

The equation (01) gives a eondition whieh ensures that the linearization of 
the equilibrium equation Q has a niee strueture whieh allows to implement 
a Newton method with tame estimates (namely, that it faetorizes into two 
first order differenee equations). See Seetion lA^ for a diseussion of dH) as a 
tool for solving Q. 

The equation dH) is hard to solve exaetly, but we will be able to develop 
a quasi-Newton method for (01). 

The main idea of the paper is that, even if we do not know how to earry 
out a KAM theory for the equilibrium equation (O alone, we ean earry out a 
KAM proeedure for the pair of equations d2l) and dll)- As diseussed in See- 
tion 11.1.31 similar things (a funetional equation supplemented by another 
auxiliary one that makes the assoeiated linearized equation solvable) have 
happened in elassieal problems in KAM theory IIMos67irEli88ll . 
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Remark 4. Notice that both the equilibrium equation (O and the factoriza¬ 
tion equations (@1) are coupled because v is an unknown in both equations. 

Nevertheless, the effect of the factorization equation on the equilibrium 
equation appears only through the counterterm cr and is a mild effect (linear 
in cr). (See Section 1431 for a detailed treatment.) On the other hand, the 
effect of the variables of the equilibrium equation is very strong. Hence, 
one can think of the pair Q and (0]) almost as an upper triangular system 
of equations. Indeed, the perturbative treatment of the pair Q and (@1) has 
a skew product structure. See Section \5J] 

1.1.3. Comparison of the method in this paper with the application of re- 
ducibility. The method of adding an extra equation so that the linearization 
of the equilibrium equation is solvable, has already appeared in KAM the¬ 
ory. 

In the theory of perturbation of lower dimensional elliptic tori, the classi¬ 
cal treatment is to try to reduce the linearized equation to an equation with 
constant coefficients [|Mos671 |Mel631 . This requires extra non-resonance 
conditions and, in principle extra parameters. See, in particular IICW99[ 
|HLY06irTre91L which study the problem of breakdown of resonant tori in 
Hamiltonian systems. 

In the present case, the situation is completely different in the details 
(since we do not seek reducibility but rather factorization into two second 
order equations) as well as in the concepts (in iMos67H . the parameters are 
related to initial conditions or the characteristic numbers of the linearized 
equation). Hence the parameter count of the present method is very dif¬ 
ferent from what one could expect from reducibility. Also the reducibility 
equations have a very different geometric meaning from the factorization 
equations. 

The factorization method has analogies in higher dimensional systems 
and in elliptic PDE IIKoz83[ IMosSSB . One can think of factorization as an 
analogue of putting the PDE in divergence form. The transformation to 
divergence form is achieved in IIKoz831[Mos88H for elliptic operators taking 
advantage of an identity (which is analogous to the Ward identities in gauge 
theory). Here, on the other hand, we have to do a KAM theorem to obtain 
an auxiliary function that gives the factorization. 

The present method gives an a-posteriori theorem and smooth depen¬ 
dence on parameters, justifies the perturbation theory and leads to efficient 
numerical implementations. 

1.2. Properties of the equilibrium equations Q and its associated fac¬ 
torization equation dH). Before embarking on the hard analysis, in this 
section, we derive some identities and symmetries of the equations which 
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are only soft analysis. This section can be skipped by readers interested 
only in the KAM methods. 

Many of the symmetries and elementary properties derived for the equi¬ 
librium equation in HdlLSZlSI lift straightforwardly to the factorization 
equation. Surprisingly, the formal perturbation theory developed here for 
the pair @ and dH) is more efficient than the perturbation theory for dU) 
alone developed in HdlLSZlSI . The perturbation theory developed in HdlLSZlSI 
was only for perturbation around integrable solutions, but the perturbation 
theory for the pair Q and dH) is developed around any solution of both 
equations ^ and dll)- See Section [51 

The expansions around zero found in jdlLSZlSl are a particular case of 
the expansions found here since the linearization of the equilibrium equa¬ 
tion around zero admits a trivial factorization. 

1.2.1. The symmetries of the equilibrium equations d2l)- We note that the 
symmetries for the equilibrium equation found in ffdlLSZlSI extend to the 
factorization equation. 

We have that if v^,, o-{q), A(q) is a solution of d2l), for any function L{q) so 
is: 

V;,(<A) = V^+LiriwM + + ^(^7)’ 

(5) o-iq) = o-{q + 

A(q) = A{q + L{q)fir,) - Lib) o-{q + i{q)Pf). 

Here we use to denote the first d-l components of {A and to denote 

the last component of (A. 

Since the symmetry dH) involves changes of arguments, giving a v,,, find¬ 
ing the L{q) that accomplishes the normalization involves solving the im¬ 
plicit equation 

(6) I{q + jAr^iiii)) + Liri) = Q 
where I{q) = Vr,(f) df. 

Applying the finite dimensional implicit function theorem, we can solve 
d^ if 7 and its derivative with respect to dq) are both small. In contrast, in 
the non-resonant case treated in nSdlL12bll , the normalization of the func¬ 
tion V could always be solved explicitly. 

As we will prove in Section jH the solutions of ^ that satisfy the nor¬ 
malization d3l) will be locally unique. 

2. Preliminaries 

To formulate the KAM results (as well as to make quantitative the Lind- 
stedt series) we need to define precisely the norms of analytic functions. 
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In this section, we collect the definitions and some standard properties of 
spaces of analytic functions. 

In Section [2m we collect several standard definitions of spaces and present 
some preliminary results on these spaces. In Sectionwe present defini¬ 
tions of the Diophantine properties we will use in this paper. In Section [231 
we present well known estimates for cohomology equations, which are the 
basis of the KAM procedure. Besides the customary constant coefficient 
equations, we study first order cohomology equations with non-constant 
coefficients in section [Z4l which were also studied in [|Her83l . 

2.1. Spaces of functions we will use. We will use a variation on the same 
spaces of analytic functions which have been used very often in KAM the¬ 
ory since [|Mos67M . We will use the same notations as in [|dlL08[|CdlL10b[ 
ISdlL12blfSdlLT2al . 

We denote by 

Dp = { ^ G |Im((A,)| < p }. 

We denote the Fourier expansion of a periodic function v(i/^) on Dp by 

v(«A) = Xi 

isZ''-! 

where • is the Euclidean scalar product in and Vk are the Fourier coef¬ 
ficients of V. 

We denote by s^p the Banach space of analytic functions on Dp which 
are real for real argument and extend continuously to Dp. We make s^^p a 
Banach space by endowing it with the supremum norm: 

llvllp = su_p Iv(l/r)|. 

ipeDp 

The spaces of analytic functions s^/p are the same spaces as in ^Mos67ll 
and that some of their elementary properties used in the argument were 
discussed in [|CdlL10allSdlL12bll . Notably: 

• Interpolation inequalities (Hadamard three circle theorem): 

(7) \\v\\ep.(i-e)p' < llvli; llvU’r^ 

• Cauchy inequalities: 

\\D’v\\p_s<C{l,d) 6-‘\\v\\p, 

• The regularity of the composition: 

Proposition 1. Let f be an analytic function in a domain ^ c C, v G 
s^p. Assume v(T^) c dist(v(Tp), C - f > 0. Then, 
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(1) / o V G £/p; 

(2) Tfllvlip < f, we have 

||/(v + v) - /(v) - /'(v)v||p < C||v|g. 

The cohomology equations we will have to consider are different from 
those studied before and we will present the results in Section [2!^ 

Remark 5. The method of proof works also for spaces of functions with fi¬ 
nite differentiability. Indeed, KAM theory is often formulated as an abstract 
implicit function theorem for the functional S acting on spaces of functions 
that satisfy some mild properties HSchbOl lZeh75ll . (The paper HCdlLlOall 
presents an implicit function theorem well suited for the method in this pa¬ 
per. ) 

In particular, the method in this paper works as well when v is considered 
in Sobolev spaces of high enough regularity. For simplicity, we will not 
formulate the finite differentiable version of the results. 

2.2. Diophantine condition. We will require that kl satisfies the Diophan- 
tine condition in : 

(8) \k-Q-m\> K\kr V /t G - {0}, m G Z. 

Here k, t are positive numbers. 

We denote by 3!{k,t) the set of kl which satisfy ([S]). We also denote 
^(t) = U^>o^(/c,t). 

It is well known that the set ^(t) with r > 1 is of full d- 1 dimensional 

Lebesgue measure in R^~^ 

2.3. Cohomology equations. It is standard in KAM theory to solve for v 
given (f) with zero average in such a way that: 

(9) v(i/r + Q) - v{f) = fif), 
where kl g S>{k, t). 

We will use equations similar to dH) frequently with some fixed frequency 
kl. To simplify our notations, we will denote v{f + kl) and v{f/ - kl) as v+(f/) 
and V-(if/), respectively. Similar notations will be used for other functions. 
We also use T to represent the translation operators, i.e. [rf 2 v]((/f) = v(if/ + 
kl). 

Estimates for ® for Diophantine frequencies were proved in [|Rus75[ 
lRus76l . The crucial point of these estimates is that the solution is bounded 
in smaller domains and that there are bounds on the solution in the smaller 
domains. These estimates have large constants if the loss of domain is small, 
but the constants can be chosen to be a power of the domain loss (tame 
estimates). See (fT^ . 
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Tame estimates (and henee Diophantine eonditions) are used in the eon- 
vergenee proofs and in the KAM theory developed here. 

Lemma 1. Let ^ e be such that 


( 10 ) 


f 


(p{i//)chf/ = 0. 


Assume that Q e 3>{k, t). 

Then, there exists a unique solution v of ^ which satisfies 


( 11 ) 

The solution v is in 

( 12 ) 


r 


v(f)df = 0. 


' for any 0 < p' < p, and 


< Cid, t) K ' (p - p') 


lip- 


Furthermore, any distribution solution of dH) differs from the solution 
claimed before by a constant. 

Iff is such that it takes real values for real arguments, so does v. 

Similarly if we consider analytic functions f e ^(T'^) satisfying fif, q)dil/ 
0, then for each q, we can solve 


(13) 


v(iff + Q.,q)- v(f, q) = 4>{f, q). 


The solution v e is analytic in {if, q) and we have 


(14) 


< C{d, t) K ^ ip-p') 


lip- 


We note that, as it is well known that obtaining v solving ® for given f 
is very explicit in terms of Fourier coefficients. If 








then, V is given by 








The above formula for the solution makes it clear that if f is real for 
real values of its arguments, so is v. Note also that if the function f is 
discretized in terms of N Fourier coefficients, the computation of the Fourier 
coefficients of v takes only N operations. It also makes it clear how to obtain 
the tame estimates (with a worse exponent) for Diophantine frequencies or 
the results for i 
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2.4. Cohomology equations with non-constant coefficients. In this sec¬ 
tion we consider a generalization of the above theory for equations with 
non-constant coefficients. These equations have also been called twisted 
cohomology equations in HHerSSH . Our strategy is similar to iHer83L but 
we will need more details about the dependence of the solutions on the co¬ 
efficient. 

We will be considering equations of the form 
(15) -I- O) - b{ij/)v{il/) = T -I- (p{ij/) 

where if/ e and a, h are fixed functions in s/p. 

We consider that (p e is the known data of the equation and T € R, 
V e ^p-s are the unknowns we seek for. We will refer to (fTSl) as “twisted 
cohomology equations”. 

We will show that, given some non-degeneracy conditions in a and b, we 
can obtain estimates for A and v. The idea, already present in [|Her83L is 
that we can rewrite the equation (fTSl) into a constant coefficient equation. 
We will go through the procedure in details because we will need rather 
detailed estimates on the effect of a and b. 


Lemma 2. Suppose \\a - l||p < r < \, \\b - l||p < r < 1 and that Q 
satisfies the Diophantine condition ([8]). Then there exist positive real valued 
functions ya,7b ^ cmd real numbers a, b such that 


(16) 


aw) = a -——, 

yaifi) 

ybi.fi) 

ybifi + Tl) ■ 


bifi) = b 


In addition, for the twisted cohomology equation (fTSl) . there exist a unique 
solution T 6 R and v e s^fip^for (fT5l) satisfying the normalization 


(17) 



vdfi = 0 


such that 


(18) 


I ll0llpllTa(y^))+llp-i5/2 
\A\ < —p - 

jjd-i yaiyb)+dfi 

\\yaybv\\p-6 < Cd^^WiA -H fi)yaiyb)+\\p V d > 0. 


Therefore, 

(19) ||v||p _5 < C||yallpl|l/rallp-5llrfollplll/7&llp-55“''ll(^ + 0)llp 


Vd>0. 
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( 21 ) 


Proof. To show (fT^ . it suffices to observe that, taking logarithms, (fT^ is 
equivalent to: 

{\oga){ilj) = log(a) + \og{ya){ilJ + Q) - logCrJC^!^) 

(log b){il/) = log(h) + \og{jb){f) - \og{yb){.f + 

which are eohomology equations with constant coefficients. Applying LemmafU 
we get solutions ya,yb £ ■s^p-eii, a, ^ e R and estimates 

II log 7a I Ip- 5/2 < C((i, T,A:)d“'"|| log a - log d||p, 

||log 7 fo||p_ 5/2 < C(J,T,/c)d~'^||logZ7-logh||p. 

Once we solve the constant coefficient cohomology equation (l20l) for 
log a, log b, log 7 q, log yb we obtain a, b, ya, yb by taking exponentials. This 
ensures that they are positive for real values of the argument. We refer to 
a, b as the average eoeffieients of the cohomology equation. 

Onee we have the solution of (fT^ . we realize that the equation (fT5l) is 
equivalent to 

(22) d( 7 a 7 fov)+ - Hyaybv) = (d + (f>)yai7b)+ 

whieh is a cohomology equation with eonstant coefficients. Let us denote 
m = yaybV for simplieity. 

When a = b, we ean solve (l22l) using Lemma [Hand get estimates (fTSl) . 

Note if a b, the equation (1221) is easier to solve sinee no small divisors 
appear. 

• Ifd ^ h, then (1221) is equivalent to 

b A + (f) 

m+- -m = ^^ya(yb)+, 
a a 


namely, m+ = ^ya{yb)+ + \m. 
Therefore, 

^bV A + f(if/ - (n + 1)Q) 


'A (b 

n=0 ' 


yaif -(n + \)Q)yb{il/ - nO.) 


is a solution of (1221) . Henee, ||m||p < C{a,b) ||(/l + (p)ya{yb)+\\p- 
• Ifd h, then (1221) is equivalent to 

a A + (j) 

^m+ - m = -^^ya{yb)+, 
b b 


namely, m = |m+ - -^ya{yb)+- 
Therefore, 

/dd + (pif + nO) 


rt=0 


yaif + nQ)ybif + {n + l)f2) 
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is a solution of (l22l) . Hence, \\m\\p < C(a,b) ||(/l + ^)ya(7b)+\\p- 
Note that when a b we do not have any loss of domain, but the esti¬ 
mates depend on a, b. To get estimates uniformly in a, b, we need to use the 
Fourier method. □ 

The choice of the parameter A so as to achieve the normalization deserves 
some discussion. In the case a = b we see that we have to choose A in such a 
way that the right hand side of (1221) has zero average. In this case, however, 
we can choose the average of the solution of (1221) arbitrarily. 

In the case a b, given any A, the solution of (1221) will be unique. Fur¬ 
thermore, the solution will be an affine function of A and, hence, so will be 
its average. 

We see that the derivative of the average of the solution with respect to A 
is 

—<m) =-=—. 

dA a — b 

In summary, in the equal average coefficient case, the equation requires 
adjusting one parameter to be solvable, but it gives back one free parameter 
of solutions. In the different average coefficient case, we do not require any 
parameter to ensure the solutions, but the solution is unique. 

In both cases, the solutions of twisted cohomology equations require as 
many parameters as they give back. This allows us to discuss the solutions 
of the equations which factorize into twisted cohomology equations and 
they require as many parameters as they give back. 

Remark 6. It will be important for subsequent applications that the esti¬ 
mates (ITSl) are formulated in terms ofy. We could have attempted to formu¬ 
late them in terms of a, but this is not practical for subsequent applications. 
On the other hand, we have formulated the normalization condition (ITTl) in 
a way that it is independent of y. This will be important in our iterative 
scheme because y will change from step to step. 

In subsequent applications, we will be using Lemma^when the a, b are 
changing. It will be important to track the changes ofy a, yt when we change 
a, b. See estimates (l69l) . 


3. The KAM theorem 

In this section, we will state precisely the main result of this paper. 

3.1. Statement of the main result. For the following theorem, we will fix 
the parameter q and omit the subscript q. 

We denote {g) the average of some function g and g = g - {g). We also 
use the notation: 

^ = Tn + r_n - 2 + dpW{{f, q) + M + m. 
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Theorem 1. Let a g such that a • j 0, 7 g \ {0} and a> £ be 
such that coa is resonant. Let W be an analytic function defined in a domain 
Q) c jjf. Take p > 0 and 0 < s < pjl. 

Denote: 

(23) (S{v, cr. A) = v+ + v_ - 2v + Wdfi, rf) + ySv) + crv + A, 

(24) ^(v, cr, c) = (-C + 2 - d^W((if,ri) +/3v) - cr)c+- 1. 

We assume: 

(HI) Initial guesses: Let (v°, c°, cr^, 4°) be an approximate solution such 
that ||(^[v°, cr°, /l°]||p < e and cr°, < e. 

(H2) Diophantine properties: There exists a matrix B £ S L(d, Z), Q. £ 
R^“\ L £ if such that Bcoa = (12,0) + L with 

\l-Ll-n\> K\ir V I £ - {0}, n£Z. 

(H3) Non-degeneracy conditions: 

||c° - Hip < Ml < 1 , 

|(r°| < M 2 < 1 , 
ll^vo|lc 2 (^) < M 3 < 1 

||v°||p < m < 1, 

where Wp is short for W((if/,rj) + /3v\ Mi + M 2 + M 3 < 1 and m 
depends on d, r, k. Mi, s. 

(H4) Composition condition: Denote 

We assume c LA and dist(,^y 0 , dQi) > 2^ > 0, where L depends 
on d,T,K,Mi,s, ^li. 

Let p'=p-s-6<p-s. Let e be such that 

€ < e*5^\ 

where e* depends on d, r, k, M, s, |yS|i and will be specified in the proof 
Then, there exist functions v*,c* £ sfp' and numbers cr*. A* G R such that 

cL[v*,(r*,A*] = 0, 

,^[v*,cr*,c*] = 0 . 

In addition, 

||v° - vllp. < C6-^^e, 

| 4 ° - 4*1 < 

|(r° - (r*| < C6-^^e, 

||c°-c*||p, <Cd-4^e, 


(26) 
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where C depends on d, r, k, M, s, |/3|i. In the rest of the paper, we will use C 
to denote any constant depending on these parameters. 

We also have local uniqueness: Suppose (vi,cri,Ai,Ci) and (v 2 , cr 2 , A 2 , C 2 ) 
satisfying 

(27) = Av 2 ,(r 2 ,^ 2 ] =0, 

(28) ^[Vi,cri,ci] = ^[V 2 ,Cr 2 ,C 2 ] = 0 , 
and the normalization condition 

(29) I Vidifj = I V 2 (ii/f = 0. 

J'fd-l Jjd-i 

If 

(30) max{||v 2 - Vi||3p/2, |cr 2 - cr 2 |, U 2 - ^il} < e*p, 
then, we have 

(31) (Vi,Ci,Cri,ii) = {V2,C2,(T2,A2). 

Finally, we have that the solution depends on the parameter q analytically. 

The Theorem [T] is in an “a-posteriori” format. Given an approximate 
solution whieh satisfies some non-degeneraey eondition, we establish that 
there is a true solution nearby. 

For the experts in KAM theory, we note that we have two parameters s, 6 
measuring the domain loss. The parameter s measures the domain loss in 
the first step and the parameter 6 measures the domain loss in the subsequent 
steps. The first step sets up several quantities whieh affeet subsequent steps 
and are estimated in a different way from the other steps. The steps after the 
first are all very similar, and the estimates are very similar to the standard 
Nash-Moser seheme, even if the iterative step is very different. 

To produee the approximate solution, we ean use a variety of methods. 

In the ease that W is small, we ean take as approximate solution (v°, c°, cr°, f 
(0,1,0,0) (the solution eorresponding to VF = 0) and then (HI) just be- 
eomes smallness eondition on W. 

We ean use the result of a Lindstedt series as approximate solutions. 
Then, we obtain a validation of the Lindstedt proeedure, and a eonsequenee 
for the (eomplex) differentiability of the solutions and the eonvergenee of 
Lindstedt series (see Seetion[5]). See also HdlLSZlSH . 

Even if we will not diseuss it in this paper, one eould take as an approxi¬ 
mate solution the outeome of a numerieal eomputation and Theorem [T] gives 
a validation of the numerieal result (provided that we eheek a few “eondi¬ 
tion numbers”). Indeed, the proof leads to an effieient algorithm. 
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4. Proof of Theorem [U 

4.1. Outline of the proof. The proof of Theorem [His based on a quadrat- 
ieally eonvergent iterative method. The eonvergenee of the method will 
be established by a Nash-Moser argument. Here we give an outline of the 
proof. 

The linearization of the equilibrium equation is a linear seeond order dif- 
ferenee equation. We will solve this linearized equilibrium equation using 
that it faetorizes into two first order differenee equations. In our ease, we 
eannot write down the faetorization in elosed form (as in [|SdlL12b[|dlL08l ). 
So we will impose an auxiliary equation for the eoeffieients of the faetor- 
ized equation. We will call this equation the factorization equation. We will 
not need to solve exactly the factorization equation at each step, but we will 
require that there is an approximate solution of the factorization equation 
with an accuracy comparable to the accuracy of the approximate solution of 
the equilibrium equation. Then, we will derive an iterative method that im¬ 
proves both of the equilibrium equation and the factorization equation. This 
is reminiscent of the procedure to study an elliptic torus and its reducibility 
at the same time. See [|Mos67[rPos89ll . Note without the parameter cr, we 
can not find coefficients satisfying the factorization equation. 

Even with this novelty, we still have difficulty that the equilibrium equa¬ 
tion and the factorization equation are coupled. We do not know how to 
establish convergence if we solve one equation after another. Therefore 
we will derive a way to solve the two coupled linearized equations simul¬ 
taneously. After we factorize the linearized equilibrium equation, we will 
face three twisted cohomology equations (difference equations with non¬ 
constant coefficients), two for the linearized equilibrium equation and one 
for the linearized factorization equation. We will solve them using the tech¬ 
nique in Section 14^21 This completes one step of iteration. 

One small technical complication is that the first step is estimated in a 
different way from the subsequent steps. In the first step of the iteration, 
we compute some auxiliary function related to factorization (the auxiliary 
functions y related to the twisted cohomology equation) and estimate them 
from the original data. In subsequent steps, even if we compute the same 
auxiliary functions, we do not estimate them from the data but we estimate 
the change on the auxiliary functions induced by the (rather small) changes 
in the approximate solution. 

Under some minor non-degeneracy conditions, we can repeat the process 
indefinitely. Finally, we will prove that the iterative procedure converges 
with suitably chosen domains under the analytic norm. This is very standard 
in KAM theory. 
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The local uniqueness will be obtained in Section |4]8] by showing that the 
linearized equation admits unique solutions using a standard technique. See 
nCdlLlObl . 

We also note that the main part of Theorem [T] has what is called “a- 
posteriori” format in numerical analysis. We show that if there is an ap¬ 
proximate enough solution that satisfies some non-degeneracy conditions, 
then, there is a true solution nearby. It is well known that this a-posteriori 
format leads to smooth dependence on parameters, bootstrap of regularity 
and several other consequences. 

Finally, the analyticity with respect to parameters(we will use either e or 
77 as parameters) will be a corollary of the existence of the perturbative ex¬ 
pansions (see Section [5T]) . Using the a-posteriori format of Theorem [Hand 
the local uniqueness, it follows that the solution is complex differentiable 
in the parameters. Hence, it is analytic on the parameters. Furthermore, 
since the perturbative expansions are the Taylor expansions of an analytic 
function they converge. 

We note that the iterative process described above leads to an efficient 
algorithm, which we have formulated in Section |4^ 

4.2. Motivation for the iterative step. Our goal in this section is to devise 
a procedure that given an approximate solution produces another approxi¬ 
mate solution with much smaller error and not much worse non-degeneracy 
conditions. This procedure is done for the equilibrium and factorization 
equations simultaneously. 

4.2.1. The equilibrium equation. We consider the initial guess (v°, cr°, 4°) 
which solves equation Q with a small error e, where ||e||p < 6, i.e. 

(32) v° + v° - 2v° + W(((A, q) + ySv°) + (r\^ + A° = e. 

The Newton procedure for the approximate solutions of the equilibrium 
equation (13^ requires to find an update (v, &, A) satisfying 

(33) D+ -I- v_ - 2v -I- dpW^fiV + (TV° -l- cr^v + A = -e, 

where we denote dfjW = jS ■ VW for simplicity and we use Wy to indicate 
Wiii/^,q)+/3v). 

4.2.2. The factorization equation. The equation (1331) is not easy to solve 
directly. Nevertheless, we impose that it can be factorized into two first 
order difference equations with non-constant coefficients. If we accomplish 
this, we can solve (1331) using the theory of twisted cohomology equations 
developed in Section IZ4l 
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Therefore we want that the operator (recall that we denote = Tci + 
T_ci - 2 + dpWyO + cr) is factorized into two first order operators: 

(34) ^ 
where 

(35) = a{ilj)Tii - b{ifj) 

(36) -c/- = c{il/) - d{i//)T-Q, 

and a, b, c, d are the new unknowns. 

A direct calculation shows 


(37) = a(^)c+(i/^)v+ - [a(i/r)J+(i/r) + (7(i/^)c(i/r)]v + b{il/)d{il/)v-. 
Hence we want to choose a, b, c, d satisfying the following equations 

= 1 , 

(38) -[a(il/)d+(ilf) + b{i/f)c{ifr)] = -2 + d/jW^o + cr^, 

b(if/)d(if/) = 1. 


Note that the problem of factorization has always many solutions. If 
£/+£/- is a solution and g is any invertible function, then £/+£/- is also a 
solution, where 


(39) 


^ = ^+g, 


We can use this non-uniqueness to impose some extra normalization condi¬ 
tion. In this paper, we will take the normalization 

(40) bix) = 1. 

With condition (l40l) we can simplify the system (1381) to (after eliminating 
a): 

(41) (-C + 2 - dpW^ - (r°)c+ = 1. 

This is how the operator ^ comes into play. Note it is a non-linear, non¬ 
local equation which we will have to solve iteratively. By assumption (HI), 
we can solve this equation with an error / using some initial guess c°, i.e. 

(42) (-c« + 2 - dpW^ - - 1 = / 

and we have ||/||p < e. 

Then the Newton procedure for (l42l) requires solving the following equa¬ 
tion for c: 

(43) - c°c + (-c** + 2 - - (r°)c+ - - dpdpW,oclv = -/. 

Then c -l- c would be a more accurate solution. 
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We replace with in the original equilibrium equation (l3^ to 

get 

(44) £/+£/_v + + A = -e. 

Remark 7. Note that the equation (l44l) is slightly different from (l3^ . so the 
(v, &, A, c) we find is not exactly the solution of the linearized equation of 
Q. But the only difference between (l44l) and (|3^ is a term /v(c° which 
we will show is quadratic in e, f (since v will be of the order ofe). 

Remark 8. To obtain quadratic convergence, the full Newton step will have 
to improve both the equilibrium and the factorization equation at the same 
time. (A moment’s reflection shows that improving one equation and then 
the other does not lead to quadratic convergence.) We will prove our New¬ 
ton procedure using (l44l) instead of (l3^ and improving the factorization 
still converges quadratically to a solution of Q. 

In the following section, we will develop a method to solve (l44l) and (1431) 
simultaneously. In Section 14.61 we will obtain estimates for the size of the 
error. Then, we will show in Section 14.71 that they lead to an improved 
solution and that the process can be iterated and converges to a solution. 

4.3. Solving the linearized equations. Our goal in this section is to solve 
(l44l) and (l43l) simultaneously. 

The difficulty arises because the equations (l44l) and (l43l) are coupled (the 
unknown & appears in both of them). Observe that since & appears in an 
affine way, we can guess that the solutions will be affine in &. This allows 
to uncouple the equations. 

When we find the corrections (v, &, A, c), we can write down the updated 
solution of the equilibrium equation @ and the factorization equation (l42l) : 

(45) (v\cr\/l\c') = (v° + V, cr° + (T,/l° + d, c° + c), 

which can be used in the next iterative step. We will show that the improved 
solution (1451) is indeed more approximate (in a smaller domain). Then the 
procedure can be repeated indefinitely and converges to a solution if the 
errors are small enough. 

Now we write 

(46) V = A + &B, 

(47) A = G + &D, 

where A and B are functions, G and D are numbers. All these quantities 
will be determined shortly. 
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Then the equation (l4^ becomes 


(48) 

(49) 


+ G = -e, 
+ D = -v°. 


The Newton equation for factorization (|4^ after substitution (l48l) and (l49l) 
becomes 


(50) 



Hence the Newton step improving simultaneously the invariance and fac¬ 
torization equations consists in solving (I4S1) . (H91) . (I50l) . Notice that this sys¬ 
tem has an upper triangular structure, we can solve (l48l) . (I49l) for A, B, G, D 
and then solve (l50l) for c and &. Once we have &, we can determine v and 
A using (|4^ and (|47]) . 

To solve (l48l) and (l49l) . we observe that each of them can be obtained by 
solving two cohomology equations with non-constant coefficients. We can 
find the zero average solution of (1481) and (l49l) . In fact, G is determined 
so that A has zero average and D is determined so that B has zero average. 
This ensures that (D) = 0. 

We note that the choice of the constants in the solution of two consecu¬ 
tive twisted cohomology equations is a simple extension of the arguments 
developed at the end of Section [ZH 

In the case that £/+, both have different average coefficients, we see 
that the average of the solution as a function of G is the composition of 
two affine functions, hence affine and if the linear part of each of them is 
not zero, we obtain that the composition of the two affine functions has 
non-zero derivative. 

In the case that has equal average coefficient, but does not, we 
choose G so that is solvable and choose the average of its solution so 
that the solution of satisfies the normalization. 

In the case that has different average coefficient but has the same 
average coefficient, we choose G so that the solution produced by satis¬ 
fies the compatibility conditions for Then, the solution of is defined 
up to a constant which can be chosen in a unique way to satisfy the normal¬ 
ization condition. 

In the case that both have the same average coefficients-which is 

the case that appears in standard KAM theory-we choose G so that is 
solvable. Then, we use the free parameter of the solution of to ensure 
the solvability of and use the free constant in the solution of si- to adjust 
the normalization. 


Therefore, we obtain the solution of (l48l) . Similar procedure can be done 
for®. 
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Now that we have solved (1481) and (l49l) we turn to solving (l50l) . As long 
as 



(51) 


whieh ean be proved under some non-degeneraey eonditions on the initial 
guesses, equation (l50l) is a twisted eohomology equation. Provided that -c° 
and -c^ + 2 - dj}W^,o - cr° are away from zero, we ean apply Lemma to 
solve (l50l) after A and B are found. 

To make the proeedure elear, we will write the full algorithm in See- 
tion 14.41 In Section 14.51 we will carefully estimate the errors after one 
iterative step to show that the errors of the equilibrium and the factorization 
equation are both reduced quadratically after performing the corrections in 
the iterative step. Finally, we will prove the convergence rigorously in Sec¬ 
tion |3]71 

4.4. Formulation of the iterative step. 

Algorithm 1. (1) Given initial guesses v°, cr^, A^, c°, set 

(2) Calculate the errors e = (aiv^, cr°, 4°) and f = cr°, c°). 

(3) Find jao and cP satisfying 

loga° = logd° -r log(y«o)+ - logCy^o). 

Also find y^o and satisfying 


logc° = logc° -r log(yco) - log(yco)_ 


(4) Compute d+,b+,ya^,yh^ (the average coefficients and the auxiliary 
functions) for and d_, h_, ya_,yb- for 

(5a) Ifd+ 4 b+,d- 4 b-, compute 


a = 'e, fi = 'l. 


SetG = Set A = a + Gfi. 

(5b) Ifd+ = b+,d- 4 b-, choose G = -(e). Set a = + G) and 

(a) = 0. Set f I = £^C^a,j32 = -^_“^1. Set A = 

(5c) Ifd+ 4 b+,d- = b-, choose a\ = -£^f^e,a 2 = Set G = 

Set A = - Ga 2 ) and {A) = 0. 

(5d) Ifd+ = b+,d- = h_, choose G = -(e). Set a = £3ff^(e - (e)) and 
{a) = 0. Set A = and (A) = 0. 

(6) Find B and D in a similar way as we found A and G. Also set 
(B) = 0. 

(7) Find & and c by solving t l5^ . 

(8) Set V = A -I- d'B and A = G + &D. 

(9) Update u. A, cr and c and repeat the steps. 
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Remark 9. When we apply repeatedly the iterative steps to obtain esti¬ 
mates, it would be advantageous in Step 3 to use the information we have 
on ja computed in the previous steps, because it leads to better estimates. 
In the first step, we can only use Cauchy estimates on the initial data. See 
Section \4J\ 

Remark 10. The iterative method described above achieves quadratic con¬ 
vergence (see Section \4. 71) . It only entails performing algebraic operations 
among functions, composing them, taking derivatives and solving the coho¬ 
mology equations. 

If we discretize a function by the values of the point in a grid ofN points 
and (redundantly) N Fourier coefficients, we note that each of the opera¬ 
tions above requires N operations either in the grid representation or in the 
Fourier space representation. If we obtain either the Fourier or the real 
space representation for a function, we can obtain the other representation 
using the FFT algorithm that requires N \og{N) operations. 

Hence the method in Section |?] achieves quadratic convergence but no 
matrix inversion (or storage) is required. It only requires 0(N) storage and 
0(N log(A^)) operations per step. 

This has already been observed in nCdlL091 in the periodic case (both 
for short range and for long range interactions) and [ISdlL12bll for the 
quasi-periodic non-resonant case. Numerical implementations in the non¬ 
resonant periodic case were carried out in iBdlL13H . We have not imple¬ 
mented the above algorithm, but we think it would be interesting to do so. 

4.5. Estimates on the corrections. Denote = -c° + 2 - dj^WyO -cr^. In 
the first step, we argue as follows. From (1211) . it is elear that all the quantities 

IITc*'IIp-i 5/45 llrjo^llp -5/4 j llp-(5/4j llp-(5/4j llyv'o||p-(5/4 and lly^o llp-(5/4 

bounded by some eonstant whieh depends only on d, r, k, M, s. We denote 
the bound by E. 

In subsequent steps, we will just assume that we have the bounds E 
for the above auxiliary quantities. These bounds will be derived from the 
bounds in the first step. 

Apply the estimate (fT91) of twisted eohomology equations for equations 
(l48l) and (l49l) . also reeall the averages of A and B both vanish we have 

(52) \\A\\p_,i2 < CE^6-^^\\e\\p = CS-^Te\\p, 

(53) IIBIU/2 < CE ^ 6 -^^\\ v \ = Cr^^Hv^Hp, 

(54) |G| < CE^Ikllp < Ce, 

(55) |D| < CE%\ < C. 

Choose m small, we ean make ||B||p- 5/2 < 1- Then 

(56) ||(-c° - d^dfsW„oclB)-X_,/, < (1 - Mi)-i(l - M 2 )-'. 
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This guaranteesL^_i ( - c° - )dif/ 0. Therefore the eohomol- 

ogy equation (l5CT is solvable. 

The estimate for (l50l) gives 

1^1 < C(||A|U/2 + ll/llp) 

(57) < + Ce 

< C6-^^e, 


and 

(58) 


llclU < CE^5-\\\A%_si2 + ll/llp) 

< C6-^^e. 


From (1^ . we have the estimate 


(59) 


l|('llp-5/2 ^ ll^llp-d/2 + l^lll■6|lp-5/2 

< C6-^^e. 


Also, from the estimates for G, D and cr, we have 
(60) |T| < C6 + C6-^^6 < Cd-^^e. 


4.6. Estimates on the improved error. Next we estimate the new error 
for both equilibrium and faetorization equation after one iteration step. 

We need the Assumption (H4) to make sure the eompositions ean be done 
in appropriate domains. Reeall that ||VFvo||^ < M. Also note v is relatively 
smaller than v*’. From (l60l) . we see that as long as e is small enough the 
range of v + D is inside the domain of W. So we ean ehoose proper ^ sueh 
that dist(e^v 0 , > / > 0. 

For the subsequent iterative step, we note that if is small enough, 
||D||p _5 will also be small. So that under smallness eondition on we 
ean ensure that v + D is well inside the domain of definition of W. The linear 
estimates are valid for all d’s, but in order to ensure that we ean apply the 
non-linear estimates, we have to ehoose d’s (the domain loss) in sueh a way 
that 6 is small. 

It is standard in KAM theory (and we will do it later in seetion 14.71) that 
one ean ehoose domain losses in sueh a way that the eomposition eondi¬ 
tion is met, sueh that S„’s go to zero fast enough, so there is still a domain 
left. Therefore, the eomposition will remain in the proper domain for all 
iterative steps and the proeedure will eonverge in a non-trivial domain. 
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Using the Taylor expansions (see Proposition [I]) and the equations (1^ 
and (021) for the initial guesses, we have: 

(61) 

S{v^+v, cr° + (T, / 1 ° + d] 

=S\}P, cr*’, /1°] + v+ + D_ - 2v + (Tv° + cr°v + &v + A - VP^o + 

=e + £/+£/_v - fv - {dpW^ + cr°)v + + cr°D + on) + d 

- VPvO + VPvO+v 

= - /v + d-v + ( 1 T, 0^0 - - d^lT„ov), 

and 

^{v^+v, cr*’ + O’, c° + c] 

=^\}P, cr°, c°] + [-(c° + c) + 2 - - (cr° + (t)]c+ 

(62) + [-C - + dpW,, - d-]c° 

=[-c - - dpW,,) - (T]C+ 

- - dpW,.)cl - dpdpW,,cln. 

Take 0 < 6 < p, we have 

(63) lK[v° + v,(r° + (T,d° + d]|U 

<\\fv\\p-s + ll^iu + IIW^^'O+v - 

<ll/llp-.5llv|U + ll<^IUI|flU + C||v|g_,. 

Also, 

(64) \mv^ + v,(T^ + &,c^+mp-5 

<(||c|U + CllvlU + ||^IU)||c|U + C\\nlJ\c\-e. 

Now it is elear the new errors are quadratie in /, v, d", c, multiplied by the 
domain loss to a negative power. This is what are called “tame estimates” 
in KAM theory. 

The final estimates for both equilibrium and factorization equations on 
the updated solutions are 

(65) IK[v° + V, (r° + d-, d° + d]||p_<5 < 
and 

(66) \\^[v^ + v,(T^ + &, c° + c]||p_5 < 

Therefore we get updates for c, v, cr and d which reduce the error of both 
the equilibrium equation and the factorization equation quadratically under 
the condition that all the bounds we derived before still hold for all iterative 
steps. We can prove the convergence following standard procedure, which 
we give in the following. 
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4.7. Proof of convergence. The main effect of the iterative step is to re¬ 
duce the error in (1441) and (1431) . But it could also deteriorate the constants in 
the non-degeneracy assumptions when it improves the solution. 

Therefore the first goal of this section is to show that the constants in 
non-degeneracy assumptions of the solution do not deteriorate much and 
that the deterioration can be estimated by the error. In this case, as we will 
see, the convergence can be proved by choosing suitable domains (this is 
the choice we will do first) for each iterative step as is standard in KAM 
theory. 

We will use subscript n to denote the quantities p, 6 and e after application 
of the iterative step n times, while we use superscript n for v, c, cr, d and B. 
We take 

(67) po = p,pi = p - 5 - do, = do • 2“" and p„+i = p„ - d„. 

Denote 6„ = ||(f(v",d”,/l”)||p,, and = ||,^(v", cr”, c")||p„. We will prove 
the following holds for all iterative steps by induction. 

(Bl): All the quantities ||v”||p„, l|c” - l||p„, \cr% ||Wp||p„, ||B"||p,., \\ycA\p„, 

llrv"IU, llry-llU, and ||(-c°-d/jd^Wpc°B)-^i|p„ are still bounded 
uniformly in n. To be precise, for any quantity A° < E, we will 
prove |A" - A”“^| < E • 2“”), therefore. A” < IE for any n. 

(B2): Denote = {(i/r, 77 ) - 1 -/?v”, {if/, rf) 6 T^J. Then dist(,^v«, d^) is 
also bounded by 

(B3): < (C6o)2""', < {C~eoT'. 

Note that as a consequence of (B3), and the choices of d„, we obtain that 
the composition assumption holds if 6q is small enough. 

The first step is already shown in Section 1431 Now we assume the first n 
steps are proved. 

We first prove (Bl). We only prove the bound 

( 68 ) llrc.«-rc"-*IU-4/4<£-2-". 

The proof of the other bounds is similar (up to changing the symbols). 

From (I 2 TI) . we have 

(69) II log Jc-r - log IU-4„/4 < Cd^ll log c" - log ||p„. 

From estimate (l58l) for c, we have 

(70) ||c”||p„ < Cd„-'^6„. 

As a result, we have 

(71) II log c" - logc"^^||p„ < C6~^^en. 

Combining (l69l) and (1741) . we have 

(72) II log Jen - log ||p„-5,./4 < Cd;^^e„. 
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Therefore we have (using that 7 c" are uniformly bounded) 
(73) || 7 c" - y^-i\\p,-sjA < C6~‘^^en. 


Note for (B2), we only need bounds for v", whieh is easy to prove. 

Now recall 6„ < {Ceof". Therefore if we choose 6 o such that 6~^'^€o is 
small enough, we can guarantee (1^ (i.e. (B2)). (B3) can be shown as 
follows 


(74) 


< (Cdo'"2'^6o)"" 
= (C2'^6o)"". 


Similar estimates also hold for €„. 

4.8. Proof of local uniqueness. Suppose we have two solutions (vi, cti, Ti, ci) 
and (v 2 , cr 2 , A 2 , C 2 ) of both the equilibrium and the factorization equations, 
which also satisfy the non-degeneracy conditions. We have 

(75) <r[vi,cri,Ti] = Av2,(T2,T2] = 0. 

Then we can write 

(76) 

^[v 2 , cr 2 , A 2 ] = (f[vi,cr 1 , Ti] -I- D^[vi, crj, Ti] • (v 2 — Vi, cr 2 — crj, /I 2 — Ti) -I- 
= DS’[Vi,(Ti,Ax\ • (V2 - Vi,cr2 - (Ti,A 2 - Ai) + = 0, 

where R^ is the Taylor remainder for S. 

Another way to read this identity (1761) is to say that Vi -Vi, 42 - 41 , 0 - 2 -cti 
is a solution of the equation 

(77) ^/+^/_(v2 - vi) -I- ( 0-2 - 0-1 )v^ -I- (T 2 - Ti) = -R^. 

Similarly, we have 

(78) 

-C+(C2-Ci)-l-(-C^-|-2-(9y3VTvi-0-')(C2-Ci)+-C+(0-2-0-l)-(9/35/3WvicJ.(V2-Vi) = -R^. 

We have shown in Section 1431 that the solutions of (177]) and (TTSl) satisfying 
the normalization condition are unique and satisfy the bounds (1571) and (l59l) . 

So we have 

(79) 

l|V2-Vi||p-P||C2-Ci||p-P|0-2-Cri| < Cp ‘^^\\R\ < Cp '‘3I|V2-Vi ||p-P||C2-Ci ||p-P|0-2-0-i |)^. 
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Since ||v 2 - vi|L ||c 2 - cilL and |cr 2 - cri| are all small, we have 


(80) 


||V2 - Vi||p <Cp 


-4r 


||V2 - Vil 


<Cp '‘’■||V2 


Vi p 

IM2 


l|V2 - 


| 3 £. 

2 


The last inequality uses the interpolation inequality (|7]). 

Therefore, as long as p“'^^||v 2 - viHv is suffieiently small (depending on 

2 

the properties of the auxiliary funetion c in the faetorization equation), we 
have ||vi - V 2 l|p = 0, whieh also implies cr 2 = cti and ^2 = Ti. 


5. Consequence of Theorem [Hand its proof 

5.1. Perturbative series around any solution. In this section we present 
methods to eompute formal power series expansions. In Seotion[5]2]we will 
show that, under Diophantine eonditions on the frequeney, they eonverge on 
a suffieiently small domain. 

Assume: 

AO: We are given a family of interaetion potentials indexed by an 
external parameter p whieh ean be eomplex. We assume that is 
analytie in both its arguments if/, rj and the parameter p. 

Al: For some parameter po we have a solution v°, cr°, d® of ([2]). We 
assume that e for some p > 0. 

A2: The operator .ifp, the linearization of the equilibrium equation 
faetorizes at po into two first order operators. 

Our goal is to find a perturbative expansion of the solutions of © in 
formal power series of p - po. Later, we will show that these perturba¬ 
tive expansions are eonvergent following an argument of [|Mos67ll whieh is 
made mueh easier by the a-posteriori format of Theorem[T] See Seetion lS© 
We seek 


(81) 


v" = v° + ^(p-po)V 

/i>0 

a^ = (r° + ^(p-po)V" 

rt >0 

d'^ = / + _^(p-po)”d" 


n>0 


in sueh a way that, when we substitute it in @ and expand (formally) in 
powers of (p - po)" we obtain that the eoeffieients of same powers mateh. 

Note that this generalizes the standard Lindstedt series, whieh is a partie- 
ular ease of the expansion in the ease that the family is just = pW and 
that we expand near po = 0. 
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If we substitute (IMl) in (l2l), expand in powers of yu - jUo, the eoeffieient of 
order n has the form: 

(82) 

- Q) + v:;(<A)(-2 + n)) + crV + + d" = 


where is a polynomial expression in v°,..., . 

The main observation is that the equation (l8^ is preeisely the Quasi- 
Newton equations dS]) whieh ean be solved by faetorization. Note that we 
get the perturbative series to all orders only assume that faetorizes at 
yu = yUo and the expansion series we get has cr = 0 to all orders. 

We note that to solve the eohomology equations without any quantitative 
estimates, as shown in fldlLSZlSL it suffiees to assume 


(83) 


lim — sup 

V-»oo A/' |i|<yv,meZ 


\n\k ■ Q, - m\ 


= 0 , 


whieh is weaker than Diophantine and, indeed weaker than Brjuno eondi- 
tions. 

We have, therefore established the following: 


Lemma 3. Under the assumptions AO, A1 and that Q. satisfies the quanti¬ 
tative conditions. 

Then, we can find formal power series as in (1^ such that for any A e N 
and any p', Q < p' < p, we have 


S’,, 


^ v"(yU - yUo)”, A”(p - PoT 


-n<N 


n<N 


< Cn,p'\p-Po\ 


N+l 


P' 


Furthermore, it is possible to find a formal power series that satisfies the 
normalization 

(vp = 0 n = l,...N. 

We refer to it as the “normalized perturbative expansion”. 

This normalized perturbative expansion is unique. 


Similarly, we ean obtain existenee of perturbation expansions for both 
the equilibrium and the factorization equations. In analogy with (1^ . we 
seek expansions 


v" = + ^(yU - yUo)"v” = / + ^(yU - yUo)”d” 


(84) 


71>0 


n>0 


= C° -p ^(yU - yUo)"c” cr^ = cr° -P - dofo-" 

n>0 n>0 

in such a way that the equilibrium and factorization equations are solved. 
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The equations for order n are 

+ v” + (-2 + + (r°)v” + vV" + A’^=Rn 

(85) 

- cy + (-c° + 2 - - cr°)cl - cV - dfidf^W'^^cy = K. 

From the proof of the KAM theorem, it’s elear these ean be solved to all 
orders assuming that faetorizes at /r = //q- The series we get do not 
have cr = 0, whieh is different from the series when we do not require the 
faetorization for all orders. 

We have, therefore established the following: 


Lemma 4. Under the assumptions AO, Al, A2 and that Q. satisfies the 
quantitative conditions. 

Then, we can find formal power series as in (l84l) such that for any A e N 
and any p', 0 < p' < p, we have 




and 




^ v”(// - poT, ^ <t"(p - PoT, ^ d"0u - pqT 


-n<N 


n<N 


n<N 


^ CN,p’\d - /2ol 


V+1 


P' 


N+l 


^ v"0u - PoT, ^ crV - PoT, ^ c”0u - PoT < Cnp \p - Po\ 

n<N n<N n<N J 

Furthermore, it is possible to find a formal power series that satisfies the 
normalization 

<vp = 0 n=l,...N. 

We refer to it as the “normalized perturbative expansion”. 

This normalized perturbative expansion is unique. 


The above results eould be interpreted in a more geometrie way. We 
have proved eonvergenee for the perturbative expansion involving also the 
eounterterms that ensure the faetorization. Given a potential W small, the 
main result shows that we ean find a loeally unique eounterterm cr, whieh 
is a funetional on W so that the equation faetorizes. The set {VF(v) + cr(W)v} 
ean be interpreted as a eodimension 1 manifold (the faetorization manifold) 
in the spaee of potentials. 

The perturbative expansions in Lemma[2]can be interpreted as the pertur¬ 
bative expansions for a path of potentials indexed by e but requiring that W 
stays on the faetorization manifold. 

In this geometrie interpretation, the eondition (ISTl) that among cr we ean 
aehieve the faetorization ean be interpreted geometrieally as saying that the 
direetion (in the spaee of potentials) given by crv are transversal to the fae¬ 
torization manifold F. Note that this faetorization manifold may depend on 

(jj. 
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5 .2. Convergence of Lindstedt series for the equilibrium and factoriza¬ 
tion equations. 

Theorem 2. Assume that the conditions of Lemma 0 hold and that Q. is 
Diophantine. Then, the normalized formal series obtained in Lemma 0 is 
convergent. 

Note that a particular case of the above result is the convergence of the 
Lindstedt series when the frequency is Diophantine. 

We need the frequency to be Diophantine because, as we will see, the 
proof uses repeatedly the KAM theorem. It seems quite possible that for 
the frequencies that satisfy ([83]) but not (|8]), for many perturbations, it is 
possible to obtain perturbative expansions to all orders, which nevertheless 
do not converge. 

Proof Without loss of generality, we assume that satisfies ([3]). 

Using the KAM theorem (Theorem]!]), we obtain that for all /i in a small 
ball centered on /io there exists a (unique) normalized solution. (It suffices 
to note that v** is a sufficiently approximate solution for all p close to pod 

Now, given any p in this small ball, we can obtain a perturbative expan¬ 
sion in powers of p - p. 

We consider the first term of the expansion vl. We remark that it will be 
uniformly bounded in || ■ ||p/. 

We note that because 

Il4[vj] -P (p -p)vl,Al -P ip -p)Ald\\p' < C\p-pf 

we can apply the KAM theorem (note that the non-degeneracy conditions 
of the KAM theorem are satisfied with uniform bounds when the ball is 
considered small enough) and obtain that there is a normalized solution 
V*, T* of (]2l) for any value of /i in a ball around p and that it satisfies 

(oO) ^ ^ 

\\A;-Al + (p-p)AX,<C\p-p\f 

Using the uniqueness obtained in the KAM we obtain that v* = v°, d* = 
dp. Hence, dSH) means that vl is the derivative at p = pof the mapping that 
to p associates v®, d® if we give v the topology in 

We recall that the Cauchy-Goursat theorem shows that any complex func¬ 
tion which is differentiable at every point, is analytic llAhl781 . This argu¬ 
ment also works for functions taking values in Banach spaces. Alterna¬ 
tively, it is not difficult to show that the mapping p is continuous (a 
quick way is to show that the graph of the map is closed because of the 
uniqueness and that, since it uniformly bounded, it is compact by Montel’s 
theorem [|Ahl78l ). 
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Once we have that the funetion is analytie, we know its Taylor series 
eonverges, but the Taylor series has to be the one given by the formal series 
expansion. □ 

Notiee as a eorollary of the dependenee on parameters we ean obtain that 
the solution is analytie in the parameter rj. 

Sinee in the physieal applieations 77 G TMs important to diseuss the peri- 
odieity in 77 of the solutions thus obtained, we remark that, if we start with 
an approximate solution whieh is periodie in 77 , we will obtain a solution 
whieh is also periodie in 77 . 

This ean be seen in two ways. One ean observe that, applying Theorem[T] 
we ean obtain solutions in small enough intervals of 77 . They will be analytie 
in these small intervals and therefore they give a global analytie solution. 
Furthermore, we observe that if the approximate solution eorresponding to 
77 = 1 is elose to the solution eorresponding to 77 = 0 , they have to agree 
beeause of uniqueness. Henee, we obtain that the solution is periodie. 

We eould also argue that, the proof is based on an iterative step and that, 
by examining the proof, all the steps preserve the periodieity in 77 of the 
approximate solutions. Henee, the KAM proeedure that we deseribe for a 
fixed 77 lifts to a proeedure for periodie funetions of 77 . 

The lifting of the problem to a spaee of funetions of 77 also gives a direet 
proof of smooth dependenee on parameters. 
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